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Abstract. We investigate the sensitivity of the solar model to changes in the nuclear 
reaction screening factors. We show that the sound speed profile as determined by he- 
lioseismology certainly rules out changes in the screening factors exceeding more than 
10%. A slightly improved solar model could be obtained by enhancing screening by about 
5% over the Salpeter value. We also discuss how envelope properties of the Sun depend 
on screening, too. We conclude that the solar model can be used to help settling the 
on-going dispute about the "correct" screening factors. 
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1. Introduction 



The solar interior is, thanks to the tremendous progress of helioseismology, known to such 
great accuracy that the Sun can in fact be used as a laboratory for physics. Examples are 
investigations about the equation of state (Dappen & Nayfonov 2000) or axion properties 
( 3chlattl et al. 1999 ). In this paper, we will apply it to the Coulomb screening of the 
nuclear reactions. 

The plasma correction to nuclear reaction rates, also known as screening, is one of 
the ingredients for stellar and solar model calculations, which has repeatedly been redis- 
cussed in the literature. The standard derivation by Salpeter (1954) discusses electrostatic 
screening in the Debye-approximation, where the electrostatic potential around an ion of 
charge Zie is 



$ = {Zie/r)exY>{-r/D), 

where D is the Debye radius {D = kT/Ane^rie 



(1) 



Ue'. electron density), within which (r <^ 
D) the electrostatic potential around Zi is reduced to 

$ = Zie/r - Zie/D. (2) 
This reduced Coulomb potential gives rise to an increased reaction rate by a factor 
'Z^Z^e^' 



f = exp 



kTD 



1 + /s = 1 + O.lSSZiZsC^P^Tg 



and 



(3) 

e = Ez^^^^ the 



Z2e denoting the partner charge in the nuclear reaction 
effective charge. 

E g. 1^ ha s been challenged in various papers (e.g. Carraro et al. 1988| ; |Shaviv fc Shaviv 



1996, 200C 



Savchenko 1999), cither developing different pictures for the configuration of 
Z2 around Zi, or by pointing out the static character of Salpeter's derivation, asking 
for a dynamical one, since we are dealing with a plasma. All these papers have been 
refut ed by other work pointing out flaws, inconsistenc ies or by r e-derivation of Eg. ^ 
(e.g. iBrueggen fc Gough 19971 [Gruzinov fc BahcaU 1998|) . Recently, [Bahcall et al. (2000a) 
have summarized five different derivations of the plasma correction to nuclear reaction 
rates at the centre of the Sun, all arriving at the Salpeter formula. At the same time 
they review several papers with deviant scr eening rates, stressing that all these papers 
arrive at different factors. In their summary, BahcaU et al. (2000a) emphasize that in the 
future proving the correctness of screening formulae should be part of the corresponding 
work. Among the work criticised was also that of one of the present authors (Tsytovich 
fe Bornatici 200C ; Tsytovich 200C ) , in which screening results in a suppression of nuclear 
reaction rates in the Sun. 

Since many of the paper in this field are not easy to understand and in particular 



their correctness not easy to be verified nor falsified, we used in Flaskamp et al. (2000) 



a different approach. Following our previous applications of the Garching Solar Model 
(GA RSOM; [Schlattl 1999) as a highly sensitive laboratory for stellar an d particle physics 
(e.g. Schlattl 199S| ; ^chlattl et al. 1999 ) we applied the screening rates of Tsytovich (2000) 



to our solar model and compared the resulting sound speed with a model with Salpeter's 
eguation. In this preliminary work we reported about changes in sound speed both at the 
solar centre and in the radiative regions around r/R^ « 0.65 {Rq is the solar radius), 
which are inconsistent with helioseismological results. 

In the present paper we are extending our analysis to general modifications of Eg. || 
and repeat it for a slightly modified derivation of the screening compared to the original 
work by Tsytovich fc Bornatici (2000)| . In Section 2 we will briefly summarize our solar 
models and describe the modifled screening values. In Section 3 we show the influence 
on sound speed, neutrino rates and global solar properties and demonstrate that, given 
the accuracy of the helioseismologically inferred sound speed, even the case of reactions 
rates without screening (/ = 1) can be excluded.]^ We then turn around our approach 

^ We note that a very similar investigation has recently and independently been published by 
Fiorcntini et al. (2000) 
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to show that a slightly improved solar model can be obtained by increasing Salpeter's 
screening by a moderate amount (w 5%). Finally, in the last section, we conclude that 
the present work demonstrates convincingly that the Sun can be used as a laboratory 
for nuclear reaction screening; this could be helpful in the future to settle the on-going 
dispute about screening. 



2. Calculations and models 

2.1. Solar models 



We cal cul ated solar models with the program and input physics as described in pchlattl 



(1999) or Schlattl et al. (1999). The physics is ve ry similar to other standard models 



such as published recently by Bahcall et al. (2000b)| . We avoid repeating details here, but 



emphasize that helium and metal diffusion are included and the present solar luminosity, 
radius and surface value of Z/ X — 0.0245 are matched ( note that Bahcall et al. 2000b use 
a slightly lower value of 0.0230; Grevesse fc Sauval 1998 ). The r esulting model paramet ers 
(Table |l|) and sound profile (Fig. |l|; solid line) agree well with Bahcall et al. (2000b) . 



For our standard model ("GARSOM") we have used Salpeter's formula. In all other 
models, screening factors were replaced as described next. Otherwise the model calcula- 
tions are identical. 



2.2. Screening factors 

Except for the Tsytovich-screening model (Fig. [^; dashed line), the screening factor / of 
all reactions of pp-chains and CNO-cycle was multiplied by a constant factor g, ranging 
from 0.9 to 1.1, such that / = (1 + /s) • g for the models shown in Fig. ^, or was replaced 
by / = 1 for the no-screening model (Fig. |l|; dash-dotted line). 

For the Tsytovich-screening model a varying screening factor depending on reaction 
and composition was used for the p + p, '^He -1-"^ He, "^He -|-^ He, ^Be -I- p, and ^Be + 
reactions. (For all other reactions standard Salpeter values were used.) These factors were 



determined by starting out with the formulae given in Tsytovich & Bornatici (2000) 
which involve the solution of multi-dimensional integrals, which depend on temperature, 
density, and abundances of the participating nuclei. Instead of evaluating these integrals 
at all points within the solar model, we did so for a selected number of mass coordinates 
for several models along a standard model evolution sequence and found that to very good 
approximation the screening factors can be expressed as a function of helium content (as 
a parameter describing both the spatial coordinate and the evolution with time within 
the energy generating core). In Eq. ||we replace ZiZ2^^ by F ■ Lij, where F is depending 
on composition and Lij on reaction ij. For the centre of the Sun we recover the screening 



factors given in Table 1 of Tsytovich (2000) when using F and Lij from Tsytovich fc 
Bornatici (2000)| . 



We note that the original results of Tsytovich fc Bornatici (2000) have been modified 



slightly by V.N. Tsytovich until the time of the computations presented here and are 
subject to further ongoing research. With respect to the original paper, the quantity 
Lij became temperature-dependent (except for e~-capture on ^Be) and some additional 
terms in expansion series had been added. Meanwhile the formalism has been improved 
to treat two-particle distribution functions, a sign-error in one term has been corrected 
for, and it appears that in addition to the dynamical screening terms the static ones, 
which correspond to the Salpeter-screening and which canceled in the original derivation 



( Tsytovich fc Bornatici 200C ) could remain, depending on whether plasma perturbations 
induced by the nuclear reactions themselves can decay in time. 

The new derivation and results will be published in a forthcoming paper (Tsytovich 
2001, in preparation). We therefore emphasize that the purpose of the present paper is 
not to investigate the Tsytovich (2000)| screening rates, because they are already under 



further development, but to demonstrate th e capability of the Sun as a laboratory for 



screening using the screening formulation by Tsytovich (2000) as an example 
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Fig. 1. Comparison of the sound speed profile for our three models using Salpeter- (soHd 
line), Tsytovich- (dashed) and no screening (dash-dotted) with the sound speed derived 
from helioseismology (see text). The shaded area and the long-dashed line refer to dif- 
ferent error estimates in the latter 

3. Results 

3.1. Sound speed 



The results of our experiments are displayed in Figs. ^ and |^ for the sound speed 
difference between models (cm) and the seismic Sun (cs). The inferred sound speed profile 
we took from Basu (1998), For the error range of Cg we show two different results: the 
extremely conservative estimation by Degl'Innoccnti et al. (1997), shown as the grey- 
shaded are a in Fig, p], and a more recent analysis about the uncertainty of the inversion 
procedure (Basu et al. 2000). In the latter paper use of SOHO-results was made and the 
authors concluded that the relative sound speed errors due to the measurements are of 
order 3 ■ 10""* (3a) and that inversion method and starting model add about equally, 
thus that a conservative error range based on this paper would be of order 1 • 10"'^ for 
0.2 < r/i?0 < 0.8. This is indicated by the long-dashed line in both figures. 

It is immediately clear from Fig. ^that the model with Tsytovich-screening is stronly 
discrepant with the seismic model and even the no-screening model is outside the conser- 
vative error range for small radii and below the convective envelope. The deviation from 
the standard model is of order the conservative error range or larger, such that it can be 



excluded. This is in particular true for the more up-to-date error estimate by Basu et al. 



(2000) in the radiative interior. 



It is interesting to note that going from Salpeter- to Tsytovich-screening there are 
systematic changes in both the core and sub-convective regions, which are anti-correlated. 
The model's sound speed is increasing in the core due to the higher central temperatures 
necessary to provide the same number of pp-reactions per second (the solar luminosity 
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constraint). At the same time Cm is decreasing in that part of the radiative zone outside 
the energy-generating core (extending to r/i?Q « 0.2) due to a shghtly more expanded 
structure. The Sun appears to be more centrally concentrated than it is for standard 
screening. Recognizing this it is only a small step to ask if the two extrema in Fig. Q 
along the standard-model can be reduced simultaneously by increasing the screening 
above the Salpeter value. 

The result of the corresponding experiment is shown in Fig. ^ for different factors 
multiplied to Eq. ^ for all reactions. An increase of the screening can indeed improve 
the situation at the discrepant bump around r/i?Q « 0.65, but leads to a comparatively 
large change in the core. An increase over Salpeter's screening by more than 10% will 
become problematic even in the case of the conservative error estimates; reductions by 
5% or more can clearly be excluded. There are two lessons to be drawn from this: 

1. agreement with the seismic sound speed profile is obtained from screening factors not 
deviating by more than -|-10/-5% from Salpeter's formula; 

2. the discrepant region below the convective envelope might not necessarily be due to 
missing physics only in that region; it should not be forgotten that the Sun is a system 
coupling different regions and that changes in the core also affect the rest of the Sun. 

We close this section by reminding that the structural changes just discussed depend 
to first order only on the screening of the pp-reaction (as we have verified in test calcula- 
tions), which is the one directly connected to the solar luminosity. All statements made 
above therefore restrict the screening of only this reaction. We also remind the reader 



that the reaction rate its elf is influencing the mod el as shown by Antia fc Chitrc (1999) 



we are using the rate by Adelberger et al. (1998), where 5(0) = 4.00 • 10"^^ MeVb. This 



value, together with that of the central metallicity {Zc — 0.021) is marginally consistent 
with helioseismic data, as shown in Antia & Chitre (1999), Figure 2. 
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Table 1. Properties of the three solar models of Fig. [T[ In the upper part neutrino 
fluxes (in s~^ • cm^^) as predicted from the models are given. In the central part the 
corresponding expected measurements in the three experiments are listed; the measured 
values are given in the first column. For the CI- and Ga-experiments the unit is in 
SNU, for Super-Kamiokande it is in units of the standard model prediction. The lower 
part shows some global properties of the initial and present Sun. Note that our mixing- 
length parameter is different from those of other standard models because of the different 
treatment of convection and atmospheres in our models (see ^chlattl 19991 ). 





Salpeter 


no screening 


Tsytovich 


pp 


5.93 • 10"' 


5.96 • 10'" 


5.96 ■ 10'" 


pep 


1.390 • 10* 


1.427 • 10* 


1.463 ■ 10* 


hep 


2.076 • 10^ 


2.107 ■ 10^ 


2.456 • 10^ 


7Be 


4.812 • lO'' 


4.677 ■ 10» 


4.023 • 10^ 


8B 


5.054 • 10^ 


5.353 • 10*^ 


5.225 • 10'^ 


13N 


5.849 • 10* 


5.004 ■ 10* 


8.101 • 10* 


150 


5.072 • 10* 


4.206 ■ 10* 


7.323 ■ 10* 


17F 


6.251 • 10*^ 


5.017 ■ 10** 


9.200 ■ 10*^ 


Sum 


6.54 • 10"' 


6.54 • 10'" 


6.54 ■ 10"' 


CI (2.56) 


7.579 


7.819 


7.790 


Ga (72.5) 


128.40 


127.05 


127.61 


SK (47.8%) 


100% 


105.9% 


103.4% 


Q 


0.9750 


0.9281 


0.8731 




0.2747 


0.2734 


0.2721 


Zi 


0.0199 


0.0200 


0.0201 




0.2448 


0.2428 


0.2403 




0.0181 


0.0181 


0.0182 


Rbcz/Re 


0.7135 


0.7163 


0.7188 




1.5707 ■ lO'' 


1.5792 ■ 10^ 


1.5947 ■ 10'^ 



3. 2. Other properties 

The Sun does not only show the influence of screening in its sound speed profile, but also 
in some other properties (Table |l|; lower part). Related to results based on helioseismology 
are the depth of the convective envelope and the present surface helium content, for 



which the det ermined values are R hcz/Rp) = 0.713 ±0.001 ( Basu &: Antia 1997 ) and Yg 



0.249 ±0.003 ( Basu &: Antia 1995 ). While the Salpeter-model nicely fits, the no-screening 
case already is at the border of the error range and the Tsytovich-model disagrees with 
the value for i?bcz and is only marginally consistent with that for Ys . Note that these are 
envelope resp. surface quantities, but the change to the physics is done deep inside the 
core! 

The central temperature Tc is increasing with decreasing reaction efficiency, as ex- 
pected. This is due mainly to the pp-reaction, which is the one closely linked to the 
luminosity of the Sun and therefore reduced screening must be balanced by higher tem- 
perature. The effect is measurable because of the rather low dependence of the reac- 
tion on temperature under solar conditions. For the same luminosity reason the flux of 
pp- neutrinos remains almost unchanged (Table |l|; upper part). Temperature sensitive 
i/-emission rates should increase more {pep and hep) , unless the influence of the reduced 
screening is dominant. Considering the ^Be flux, the transition to no-screening reduces it 
by 2.7%, but the further step to Tsytovich- screening leads to a reduction of 16.3%. The 
^B-source, however, actually increases due to the increased temperature and its extreme 



temperature sensitivity. Note that in Flaskamp et al. (2000) we reported about a strongly 
reduced *B-flux. This result was wrong, because we erroneously had not modified the 
screening of the ^Be + e~ reaction, which does not appear explicitely in our nuclear 
network, but only as a rate relative to the ^Be + p-capture. 

As a consequence of the modified neutrino emission rates, the predictions for the 
neutrino experiments change. Most interesting is the case of Super-Kamiokande, which 
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the Tsytovich- screening model predicts to measure even more neutrinos than the standard 
model does, although the rate is much stronger suppressed (see Tsytovich 2000| , Tabel 1) 
than the pp-reaction. However, due to the luminosity constraint and increased central 
temperature, the modified screening effect is overcompensated. 



4. Conclusions 

We have demonstrated that the Sun has become a sensitive laboratory to investigate the 
correct plasma screening for its own nuclear reactions. The high precision of the seismic 
Sun does not allow for deviations of more than a few percent from the original formula by 



Salpeter (1954). A slightly enhanced screening factor improves the agreement with the 
seismic sound speed profile. We emphasized that not only the central parts, where the 
influence of changes to the nuclear reaction rates is immediate, but also the outer layers 
will have a modified structure pointing out discrepancies. We have to recall, however, that 
this is mostly due to changing the screening of the pp-reaction. Reactions in the higher 
pp-chains only will not influence the models significantly. Because of these reaction rates, 
in contrast, the neutrino fluxes will change, not leading to the once proposed solution of 
the classical neutrino problem (which, as is known by now, cannot be solved by changing 



reaction rates only; Hata et al. 1994), but instead to even larger discrepancies, as we 
found for the Super-Kamiokande prediction of the Tsytovich- screening model, although 
this modified screening suppresses the reactions at fixed temperature. 

One should keep in mind that physical pictures deriving alternative screening rates 
should in principle also be used to rederive the configurational effects in the equation of 
state and opacities to arrive at consistent physical model input. Our present work thus 
must remain incomplete in this respect, but we trust that the success of the standard solar 
model with Salpeter's screening is unlikely to be due to a conspiracy of using erroneous 
physics for three major model ingredients. 

We believe that the Sun itself will help to settle the dispute about the true amount 
of screening of nuclear reactions. It has become a truely invaluable celestial laboratory 
for physics. 
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